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1. Introduction 

The Morse potential [1] is one of the simplest examples of the Natanzon potentials [2] 
which has finite number of bound-states. It was introduced by P. M. Morse in 1929 as 
a model for describing the vibrating energy of a diatomic molecule and takes the form 



where x = (r — r )/r and a = ar . Here Vq represents the dissociation energy of 
molecule, ro is a bound length and a is a parameter to control the width of the potential. 

However, the radial Schrodinger equation of this potential for arbitrary non-zero 
Z-states (/ 7^ 0) cannot be solved exactly unless by using various approximate schemes. 
Among them, the most convenient one is the Pekeris approximation [3] where the basic 
idea of the method is so that the effective /-dependent potential keeps the same form 
as the potential with I = (s-state). Using this approximation, the rotating Morse 
potential has been solved using some useful methods and approaches including the 
supersymmetry (SUSY) [4], the variational approach [5], the Nikifovov-Uvanov method 
(NU) [6], the shifted 1/iV-expansion (SNE) [7], the modified shifted 1/iV-expansion 
(MSNE) [8], the asymptotic iteration method (AIM) [9], and exact quantization rule 



In the present paper, we sought to tackle the related rotating Morse bound-states 
problem differently from other methods approached before by performing the inverse 
contour representation (ICR) [IT] in the framework of the contour Laplace transforms, 
in order to solve the radial Schrodinger equation. The original statement for the use of 
ICR approach demanded a simply knowledge of the integral representation of the Euler 
beta function based on Cauchy's integral formula. This last restriction may be relaxed 
by a creation of a barrier, a cut line, joining the points and 1, which represent the 
finite limit integrals of the beta function. We will see that ICR approach is conceptually 
simple and can be extended to others potentials. 

This work is organized as follows. In section 2 the contour Laplace transforms of 
the radial Schrodinger equation with the rotating Morse potential is outlined within 
Pekeris's approximation. The contour of integration C is chosen so that the wave- 
function vanishes at the end points. Section 3 is devoted for obtaining the exact bound- 
states for the rotating Morse potential using the properties of ICR. In order to verify the 
effectiveness of our method, the energy spectra are calculated numerically for different 
diatomic molecules and compare with those obtained by others methods. In the last 
section, we do our final conclusion. 

2. Vanishing boundary condition on the wave-function 

The radial Schrodinger equation for the rotating Morse potential can be expressed as 



VuXx) = V (e 




(1.1) 



(EQR) pro]. 




(2.1) 



Inverse contour representation as a solution 



3 



where m is the reduced mass of the molecule and Kfr.( r ) is the effective potential given 
by 

V eS Xr) = V l (r) + V M .(r) 

= ^ + + V (e~ 2ax - 2 e~ ax ) . (2.2) 
2mr 

It is well-known that ( 12. lj) associated to the effective potential ( 12.2ft cannot be 
solved analytically for an arbitrary I ^ due to the centrifugal term. In order to obtain 
approximate analytical solution, we most deal approximately with the centrifugal term 
using the Pekeris approximation [HHHI] to transform the /-states rotating Morse potential 
to the s-state one. 

Then the centrifugal term in (12.21) can be approximated as [HI HO] 

1] f « 7 (oo + a x e— + a 2 e~ 2ax ) , (2.3) 
2mr z ' 

and expanding the right-hand side of (12. 3p around x = 0, we obtain the following 

parameters 

l{l + l)h 2 , 3 3 4 6 1 3 /0/1 . 

7 = — ^ o - ' fl o = 1 1 — o, oi = o' fl 2 = 1" ~^ 2 .4) 

2mrQ a ar a ar a a z 

Substituting (12.21) and (I2.3p into (12.11) the one-dimensional Schrodinger equation 
for the new rotating Morse potential Kff.( r )> with respect of the variable x, is given by 
d 2 



^ 2 _ p* e -"»* + /3fe- ax - e 2 ) 7^(a;) = 0, (2.5) 
where 

/3 2 = A 2 (Vb + 7« 2 ) , /9? = Ao (2^o - 7«i) , e 2 = -A 2 (E U;l - 7 a ) ,(2.6) 
with Ao = roV^m/h. Invoking the change in variable as well as in wave- function as 

i = e~ ax and K^x) = CHO, (2-7) 
and letting k = e/a allow to reduce (I2.5P into 



HO = o, (2.i 



where (12.81) accepts a regular point at £ = and an irregular one at £ = oo. 

Now let us consider the function expressed in (12. 8p as an integral of the general 
form [12] 

HA0 = f c Q P , q (t)^dt, (2.9) 

where Q p , q (t) is an unknown function to be determined. 

In ( 12. 9p . C is the contour of integration in the real axis used to impose the vanishing 
boundary condition for the eigenvalue problem of (12. 8 p and which does not depend on 
variable £. The integral ( 12.9P is often called the contour Laplace transforms. 
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Then by applying the n -derivative of J- p , q {Q with respect of £ and multiplying 
the resulting equation, i.e. (I2.10p . by £ and performing derivation by parts, we obtain 

^(0 = £t n Q M (t)e^dt, (2.10) 

C^S(e) = {t"6 P , g (t)e^} c -£ d[tn ^ (t)] e*dt, (2.11) 

where the symbol 01c denotes the increase of the function Y(t; when £ describes 
the contour C. Substituting (j22J), (PXTOl and fTZTTl) into ( ET5)I . we have 



« 2 / j < 1 



£ { { f - §) s » (t) + K 1 - !) < - §] s - (t) } e5 ' di = °- (2 - i2) 

where the terms in (12.121) must vanish individually. Indeed the contour of integration C 
must be chosen so that the first term in (I2.12p vanishes as well as the integrand at the 
end points; these considerations lead to solve the differential equation on Q p , q (t) and we 
obtain 

Q P , q {t)=K, q {^-tj (^ +t Y > ( 2 - 13 ) 

where Af p , q is the reduced normalization constant. Here the parameters p and q are not 
integer and defined positive 

e 1 p\ 2 e 1 A 2 . 

a 2 2ap 2 a 2 2ap 2 

Then inserting ( 12 . 1 3[) into (I2.9p . we obtain the solution of our problem which can 
be recast as 

Fp«{S)=Mp,A (--*)" (- + t Y ( 2 - 15 ) 



jq \ a j \a / 
and as a consequence the contour C, by virtue of (12.12p and (12 . 1 3j) . must verifies the 
condition 

HXH* 6 !^' (2 ' 16) 

Substituting the new change of variable t — ^(2 z — l)/a into (12. 151) we found that 
the integrand has a branch points at zq = and Z\ — \ and then the product z v (1 — z) q 
in ( 12.16P will vanish for zq and z\. The integrand into ( 12. 15ft is therefore single- valued 
for the contour encircling both branch points; in other words, taking the line segment 
joining z = and z\ = 1 as a cut line. 

Using the substitution t = ^(2 z — l)/«, the function J~ p , q (£) can be rewritten as 

•W0 = e-^/ P , 9 (0, (2.17) 
and ( 12 . 151) becomes 

f p , q {Z)=M p , q £z p - 1 0.-z) 9 - 1 e*'dz, (2.18) 
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with A = 2(5 2 /a. Now (I2.18P is in a suitable form for applying the inverse contour 
representation as reviewed in the next section. 



3. Inverse contour representation and bound-states solutions 

The inverse contour representation jTT] is an approach which may be used to advantage 
in handling differential equations by means of contour integration. 

To illustrate the technique, let us derive a familiar expression of elementary 
functions used. Let f(z) and g(z) be analytic functions of z on a closed contour C, 
then the last couple of functions are related to (" 12 . 18[) by an expression of the following 
form 

= ^~<f (V^V) 9 M (*) d*. (3.1) 

and 



2 7ri J c \z — £ 



i 

n-l 



g M (z)= / (l-x) n - l f; iq (zx)dx, (3.2) 
Jo 

where n is a positive integer and z n / (z — £) n is an analytic function provided with a 
cut line along the 7-axis joining the point and z. (13 .11) and (I3.2p are called direct and 
inverse transforms, respectively. 

We proceed here below to prove that these latter constitute a complete set of 
transformations; this step of evaluation can be performed consistently with the help of 
residues at the pole z = 00 (or u — 0) of order k + 1, where z = 1/u. 

Proof. We shall attempt to demonstrate this assertion by substituting in (13.11) and (13.21) 
a power series with undetermined coefficients bk and c^, i.e. 

00 00 

fpM) = J2 b ^ k and 9 P , q (o = E c ^ k - ( 3 - 3 ) 

fc=0 fc=0 

To do this inserting (13. 2p and the first summation of (I3.3P into (13. ip and 
interchanging the order of integration by integrating with respect of a variable x 
leading to express the Euler beta function, and comparing the final result with the 
first summation of H 3 . 3 1) we get an useful and interesting formula 

? = —B(k,n) & (3.4) 

2tti Jc {z ~ 

where B(k,n) = T(k)T(n)/T(k + n) is Euler's beta function [13]. The useful way to 
perform this integral is to use the calculus of residues; to this end the change of variable 
z = 1/u transforms the contour integral in (13 .4p into 

- — / (1 "I?"" d ^ (3-5) 

which have a pole at u = of order k + 1 . It is obvious that the minus sign in (13. 5p comes 
from the clockwise integration due to the last change of variable. Then the residue is 
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calculated straightforwardly, one obtains 

I u k+1 J u ^ k\ u^o du k 

= ^|V, (3.6) 

where (n)k = n[n + 1) • • • (n + k — 1) = r(n + k)/T(n) are the often-used Pochhammer 
symbol [13J. This identity proves (13. 4p and thus our assertion about functions f P) q(z) 
and g p , g {z). □ 

The reader should notice that (13. ip can be rewritten in the Cauchy integral theorem 
and identifying this latter to ( 12 . 18[) leads to the identity 

(w _l )!de -i ?Wfl = *P« f c ^ (1 - ^ ^ dz. (3.7) 

Substituting the second summation of (13.31) into (13. 7p and performing the (n— l)^ 1 - 
derivative and expanding exp [A£z] in the development of the Taylor series, we obtain 

f (n - a™ V r (p + fc + 1 ) r ^+1 
^ r H r ( fc + 2 ) 



K,9 A " +1 ( / zP+k (- 1 ~ ^ dz ) 



(3i 



r(fc + 2) 



fc=0 

where T refers to the Euler gamma function. The coefficient into the first summation 
in (13. 8 j) yields a two-term recursion relation identifying bk+i of ( 13.31) to i.e. 

_ r(p + fc + i) 

6fe+1 " A T(n)T(k + 2) Cfc ' (3 ' 9) 
while the comparison between the two summations leads to evaluate M p , q . 

Indeed since now z G [0, 1], then the integral calculation by means of Euler's beta 
function is straightforward which yields 

ffo) WTqj bk - (3 - 10) 

Substituting ( 13. 10|) into the first summation of (13. 3p and by introducing the beta 
function, B(p,q), inside the summation, we get the solution of our problem which is 
given by 



/ M (o = X>£* 



k=0 



= Af p>q B(p, q) 1 F 1 (p;p + q;\Z), (3.11) 

where \Fi{p ;p + q; \£) is the confluent hypergeometric function [13]. Note from (12. 14j) 
that p + q = Ik + 1 with k = e/a. 
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However it is well-known that if p equals zero or a negative integer, the series 
in ( 13.1 ip terminates and the confluent hypergeometric function becomes a simple 
polynomial. We limit ourselves here to the second case, i.e. p = —n. For this reason, it is 
convenient to redefine the confluent hypergeometric function in terms of the associated 
Laguerre polynomials xF x {-n ; 2k + 1 ; AO = n\ T(2k + 1)/T(n + 2k + 1) l1 2k) (A0, and 
substituting f[2TTTj) and fl3TL|) into (EZ} one gets 

K n M = C p , q Ce~^L^\XO, (3.12) 
where C PtQ = N p , q B(p, q) n\ Y(2k+1) /T(n + 2K+ 1) is the correct normalization constant 
to be determined. 

However from the first equation of (12.141) . i.e. 

a 2 2 afo 

we read off the bound-state energy levels where n is the vibrational quantum number. 

Indeed inserting the parameter e as defined in (12.61) into (I3.13p . the corresponding 
energy eigenvalues are then given by 

2 

+ 7^0- (3.14) 



En,l — — o 2 

Smr^ 



h 2 \ft 



it ~ (2n + l)a 

P2 



The next step is devoted to evaluate, with the help of (I3.12p . the correct 
normalization constant C p ^ q using the normalized condition 

+oo 

TZ* nl (x)TZ n ^(x) dx 

oo 

r+oo i 

= / K,i(0—, n n,m<iC (3-15) 
Jo <*£ 

This can be further written as 

C 2 m / <P 2K ~ l e-t> [L™{4>)\ d<P = a\ 2 \ (3.16) 
Jo 

where <fi = a £ and from which we obtain 

_ /20A"« / 2 t r(n+l) 

c -"(— J Vr(f + B + i)' (317) 

and the radial wave-function for the rotating Morse potential is given by 

In the derivation of (13. 17ft the correct normalization constant C p>q can be derived 
from the general expression 



S-OO 





A, T{n-k-v) T(fx + u + k + l) 



X ^ ^ r(-fc - i/)r(n - fc + i)r(fc + i)r(n + fc + i)' M ' m) 



with K(/i + i/ + 1) > 0. To prove it, write one Ln (£) in (I3.19P as the power series; i.e. 

/fc=0V - 1 / (n-k)\(k+n)\ r(fc+l) 



^(0 = ELo(-l) fc J:SLMr£i) and then use Eq.(7.41411) of Ref. [J5]. 
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4. Numerical calculations 

In tables 1-4, respectively, we give some numerical applications for four (4) diatomic 
molecules, e.g. the H 2 , CO, HC1 and LiH using f !3.14j) and then compare with the 
findings of SUSY, Variational, NU, SNE, MSNE, AIM and EQR. It is found that the 
results obtained by our approach agree with those cited above. The well-known values 
of relevant potential parameters are taken from Ref . [9] . 



Table 1 


. Energy eif 


;envalues (in eV) for the H2 diatomic molecule. 




n 


I 


ICR (present work) 


SUSY [4] 


AIM [9] Variational [5 


SNE [7] 








-4.47600 




-4.47601 


-4.47601 -4.4758 


-4.4749 




5 


-4.25879 




-4.25880 


-4.25880 -4.2563 


-4.2590 




10 


-3.72192 




-3.72193 


-3.72193 -3.7187 


-3.7247 


5 





-2.22051 




-2.22051 


-2.22052 


-2.2038 




■5 


-2.04353 




-2.04353 


-2.04355 


-2.0525 




10 


-1.60389 




-1.60389 


-1.60391 


-1.6526 


7 





-1.53743 




-1.53743 


-1.53744 


-1.5168 




5 


-1.37654 




-1.37654 


-1.37656 


-1.3887 




10 


-0.97579 




-0.97578 


-0.97581 


-1.0499 


Table 2 


. Energy eif 


;envalues (in eV) for the CO diatomic molecule. 




n 


I 


ICR (present work) 


EQR [10] 


Variational [5] NU [6] 


SNE [7] 








-11.0911 




-11.0915 


-11.093 -11.091 


-11.091 




5 


-11.0839 




-11.0844 


-11.084 -11.084 


-11.084 




10 


-11.0649 




-11.0653 


-11.0653 -11.065 


-11.065 


5 





-9.79475 




-9.79519 


-9.795 


-9.788 




5 


-9.78791 




-9.78835 


-9.788 


-9.782 




10 


-9.76967 




-9.77011 


-9.769 


-9.765 


7 





-9.29876 




-9.2992 


-9.299 


-9.286 




5 


-9.29204 




-9.29248 


-9.292 


-9.281 




10 


-9.27413 




-9.27457 


-9.274 


-9.265 



5. Conclusion 

We have investigated a method of inverse contour representation in the context of a 
contour Laplace transforms in order to solve the bound-states problem for the one- 
dimensional Schrodinger equation related to the rotating Morse potential. The main 
purpose of the approach developed here is to investigate how this method works. In 
general the determination of the inverse contour transform is the main problem using 
the integral transformation, unless if it is represented by a known function like here, 
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Table 3. Energy eigenvalues (in eV) for the HC1 diatomic molecule. 



n 


I 


ICR (present work) 


SUSY [4] 


AIM [9] 


Variational [5] 


MSNE [8] 








-4.43537 


-4.43556 


-4.4356 


-4.4360 


-4.4355 




5 


-4.39663 


-4.39681 


-4.3968 


-4.3971 


-4.3968 




10 


-4.29389 


-4.29408 


-4.2941 


-4.2940 


-4.2940 


5 





-2.80490 


-2.80508 


-2.8051 




-2.8046 




5 


-2.77193 


-2.77211 


-2.7721 




-2.7718 




10 


-2.68455 


-2.68473 


-2.6847 




-2.6850 


7 





-2.25686 


-2.25703 


-2.2570 




-2.2565 




5 


-2.22619 


-2.22636 


-2.2263 




-2.2262 




10 


-2.14497 


-2.14513 


-2.1451 




-2.1461 



Table 4. Energy eigenvalues (in eV) for the LiH diatomic molecule. 



n 


I 


ICR (present work) 


EQR [10] 


AIM [9] 


MSNE [8 


NU [6] 








-2.42876 


-2.42886 


-2.4289 


-2.4280 


-2.4287 




5 


-2.40123 


-2.40133 


-2.4013 


-2.4000 


-2.4012 




10 


-2.32873 


-2.32883 


-2.3288 


-2.3261 


-2.3287 


5 





-1.64763 


-1.64772 


-1.6477 


-1.6402 


-1.6476 




5 


-1.62368 


-1.62377 


-1.6238 


-1.6160 


-1.6236 




10 


-1.56065 


-1.56074 


-1.5607 


-1.5525 


-1.5606 


7 





-1.37748 


-1.37757 


-1.3776 


-1.3682 


-1.3774 




5 


-1.35496 


-1.35505 


-1.3550 


-1.3456 


-1.3549 




10 


-1.29572 


-1.29581 


-1.2958 


-1.2865 


-1.2957 



e.g. the Euler beta function. In fact using Cauchy's integral formula, the beta function 
may be used to advantage in the evaluation of the inverse contour representation which 
leads to the application of the calculus of residues. It is perhaps worth noting that this 
approach was successful and relatively easy, because there is not a specific choice for the 
contour of integration C in which a lot of the possible types of contours are unlimited. 
We have shown that the inverse contour representation opens up the way for another 
and new derivation of the bound-states for the rotating Morse potential. 

Finally it would be interesting to complete our approach for different types of 
potentials belonging to the hypergeometric differential equation, e.g. Poschl- Teller and 
others mentioned in |16j . 
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